Entanglement between a Two-level System and a Quantum Harmonic Oscillator 
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The entanglement between a Pauli-like two-level system and a quantum harmonic oscillator en- 
hanced by an interaction between them and a (5-pulse sequence is studied, with the decoherence 
due to their coupling with a Markovian bath. Without the Markovian bath, the entanglement is 
enhanced to maximum possible values. With the Markovian bath, the entanglement is enhanced up 
to some time and then dissipated, with the system in thermal equilibrium with the Markovian bath 
after a very long time. The time for achieving the maximum entanglement shows discontinuous 
jumps over the parameters of decoherence. 
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I. INTRODUCTION 

Entanglement is a very important quantum character- 
istic that has no classical counterpart. Einstein, Podol- 
sky and Rosen has proposed the famous EPR paradox 
in the context of two moving particles as a function of 
their coordinates and momenta, which are not commuta- 
tive py. Later, Bohm proposed the paradox in the con- 
text of two spins in terms of their spin in different di- 
rections 0. In these two examples of the illustration of 
entanglement, one is the entanglement between two con- 
tinuous systems and another between two discrete sys- 
tems. The entanglement between two discrete systems, 
especially two-level qubit like systems has been studied 
extensively with error-correction protocol. The entangle- 
ment between two continuous systems has been studied 
in depth in recent years Q. 

The entanglement between discrete systems and con- 
tinuous systems is more complicated and it contains 
much richer behaviors. Examples include an atom (two- 
leveled) coupled to a bosonic field [ij, or a system con- 
taining a superconducting charge qubit and a nanome- 
chanical resonator [^. The entanglement between 
two discrete systems and between two continuous sys- 
tems can be exchanged [7] with Jaynes-Cummings inter- 
action Q, and arbitrary entangled states between two 
continuous system can be generated arbitrarily with the 
help of a two-level system Q . On the other hand, a two- 
level system which are coupled to a bosonic dissipative 
environment are actually such a system. Sometimes an 
environment can be used to both enhance and dissipate 
the entanglement between another two systems to a cer- 
tain extent ;4| [IHl [HI- 

In this paper, the entanglement between a Pauli-like 
two-level system and a quantum harmonic oscillator is 
studied, while they are coupled to a common bosonic 
dissipative Markovian bath which induces decoherence. 
This will be organized as follows: in section |TT1 the model 
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of such a whole system and the corresponding master 
equation are discussed. The different approaches of the 
quantification of entanglement are introduced in section 
[inland their advantages and disadvatanges are discussed. 
In section IIVI without the presence of a Markovian bath 
and with the initial state being a Fock state |0) and a 
thermal state, the time-evolution of state is calculated 
and its entanglement is measured. In section|Vl the pres- 
ence of a Markovian bath is considered. 



II. 



MODEL 



In this system, there are a Pauli-like two-level system, 
a quantum harmonic oscillator and a thermal Markovian 
bath. Therefore the free Hamiltonian of this system reads 

k 



(1) 



The first term is the energy of the oscillator. The bosonic 
operator a (and a^) is the annihilation (and creation) 
operator for this oscillator. The second term is the energy 
of the two-level system with difference between the two 
energy levels given by e^. The third term is the energy 
of the bath with many degrees of freedom given by the 
momenta k's, with annihilation (and creation) operator 
being 6k (and 6k^)- 

The coupling between the two-level system and the 
oscillator is given by the following interaction Q : 



Hi = --^{o. + a))(Tz, 



(2) 



where Aq is the coupling parameter. Such interaction can 
be found in many systems, such as a system consisting of 
a nanomechnical resonator and a single Cooper-pair box. 

With only the interactions ([2]), the system would just 
oscillate between spin-up and spin-down state, as in Rabi 
oscillation. The entanglement would then oscillate be- 
tween the lowest and highest values. In order to enhance 
the entanglement, a (5-pulse sequence is added upon the 
two-level system as 
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T-Lt = ^^y^J(t-nTo)cra;, 



(3) 
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for n is any positive integers and tq = This 5-pulse 
sequence gives the effects of putting —ia^ in the calcula- 
tion, and it flips the spin of the two-level system. This 
delta function is an approximation of such flipping pulse 
and is valid is such a pulse is much stronger than the 
energy spacing of the quantum harmonic oscillator fkuo 

i- 

Both the two-level system and the oscillator is coupled 
to a Markovian bath by the following interaction: 

k 



With the density matrix of the whole system in the 
interaction picture given by p(t), its time evolution is 
given by the following dynamical equation 



(11) 



where giu and (72k are coupling constants. These cou- 
plings give rise to the dissipative effects and decoherence. 
The total Hamiltonian is given by 



Ti. = Ti.Q + Til + Tit + 'Hb- 
In interaction picture, the Hamiltonian is given by 



(5) 



V{t) — exp [ i-^t ] Hexp 



.^0 , 



By calculation, V{t) can be read as 

v{t)^njc{t) + nt{t) + nB{t), 



where 



-Hiit) 



Ao 



n 

HBit) = ^y^gik(&k^q'- + cr+bk) 

k 



(6) 

(7) 

(8) 
(9) 

(10) 



However, the entanglement between the two-level sys- 
tem and the oscillator is the concern, while the thermal 
bath is considered because of its decoherence effect. To 
facilitate the analysis, we consider the reduced density 
matrix consisting only of the two-level system and the 
oscillator, while all the degrees of freedom of the bath 
are traced out, written mathematically as 



Ps{t)^TrB{p{t)). 



(12) 



While tracing out, the bath term (flOl) is expanded up 
to second order. Suppose the bath is stationary and the 
dissipation does not depend on the memory, i.e., the dis- 
sipation is said to be Markovian. Markovian approxima- 
tion is taken and finally a master equation is derived . 
So the final equation is given as 



dt 



iaoiJQ [e~''^°\aa,psit) - p,it)aa,) + e""°*(aV,p,(t) - psit)a^a,)] 
T 

-na-{a-a+ps{t) - 2a+ps{t)(T^ + ps{t)a^cF+) 

r 

-{n^ + l)-{a+a^ps{t) - 2a^ps{t)a+ + ps{t)a+(T^) 

Q 

-nr-{aa) ps{t) - 2a) ps{t)a + ps{t)aa)) 

Q 

-[fir + l)-{a)aps{t) ~ 2aps{t)a) + ps{t)a)a) 

Y.^^^- ^^o) + e-'**a+) ps{t) - ps{t) (e***a_ + e-'**a+) 



(13) 



where 



and ao is related to the coupling by 



exp 



(14) 



Ao 
2hjjQ 



(16) 
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With the master equation, the time evolution of the 
system consisting only the two-level system and the os- 
cillator can be evaluated. 

III. QUANTIFYING ENTANGLEMENT 



well. Yet the negativity is still a good measure of entan- 
glement for mixed states. 

In this paper, both the participation ratio and nega- 
tivity will be used as a measure to measure the entangle- 
ment. 



Entanglement is a quantum property that bears no 
classical analogue. For a system in its pure state, it is 
said to be entangled if and only if there are more than 
one terms in the Schmidt decomposition of the quantum 
state where the Schmidt decomposition is given by 



E 



'A^|a„)|6„), 



(17) 



where a^'s and &„'s form two orthogonal sets of kets in 
two systems. The quantification of the entanglement is 
often given by von Neumann entropy. While the entropy 
is not advantageous for computation, the participation 
ratio [14] is often used for the sake of computation. The 
inverse participation ratio is first used as a measure of 
localization of a wavefunction [l5'| , but it is extended to 
deal with the mixedness of a density matrix. It is defined 



K = 



1 



1 



(18) 



which is the reciprocal of the purity for pr is the reduced 
density matrix for one of the subsystems. It is a reason- 
able measure because the subsystems of a pure entangled 
state is a mixed state. However, this runs into the prob- 
lem when the system itself is in mixed state, giving rise 
to the situation that the subsystem becomes mixed even 
if the system is not entangled. 

For a system in its mixed state, it is said to be entan- 
gled if and only if its density matrix can be expressed as 
a finite sum of tensor products of the subsystems: 



P - 



(19) 



One of the ways to determine whether it is entangled is 
introduced by Peres [lB|. It is proved that if the partial 
transpose of such a state consists of negative eigenvalues, 
then the state is entangled. The quantity negativity is a 
measure of the amount of these non-negative eigenvalues, 
defined mathematically as [l3] 



Af{p) 



\\P 



PT\ 



1 



(20) 



where p^"^ denotes the partial transpose of p, and \ \A\\i 
denotes the tracenorm of the operator A, i.e., ||A||i = 
Tr{V AA^). A non-entangled state have a negativity 
equal to zero. However, the converse is not true except 
for composite systems of dimension 2x2 and 2 x 3 as 
shown by the Horodecki family [l^. It often happens 
that an entangled state has a negativity equals zero as 



IV. ENTANGLEMENT WITHOUT 
DECOHERENCE 

Consider the system without the effect of decoherence, 
i.e., F = and C = 0. In that case, the entanglement be- 
tween the two-level system and the oscillator will achieve 
a maximum eventually. The state can be analytically ex- 
pressed at time t ~ utq for n being any positive integers, 
while can be numerically computed at all the other times. 

Suppose initially the state is a Fock state |0) given as 



1 



|V'(0)) = -^(| T) 



i))|0). 



(21) 



Since it is a pure state, it can be solved using the 
Schrodinger's equation 

[Hi{t) + Hmm))^^f^^^m))- m) 

And therefore, for utq < t < {n+ 1)tq, 



iV'(t)) exp 



dt'[Hiit')+Htit')]\\i;{0)) (23) 



where 



= (I T>| - {-ir&{t,nTo)]) + I i>|(-l)"a(t,nTo))) , 



a{t, nro) = 2nao + ao{l - e''"«(*-""«)), (24) 



and the coherent state is |a) = D{a)\0) where D{a) is 
given by (|A2[) . After the interaction between the two 
systems, an obvious entanglement is established between 
the two subsystems, since the level of which state the 
two-level system is occupied is correlated to the signs 
of the coherent state. To measure its entanglement, its 
participation ratio and negativity are calculated. The 
participation ratios for both subsystems are 



K 



l + cxp(-4|5(t,nro)|2) 



(25) 



And this is plotted as shown in Fig. [T] The negativity 
has to be figured out by first finding the Schmidt decom- 
position of which is 



'l + e-2|<i(t,nro)|2 

IV'W) = \ ^ — ^ l+>l/+W) 



'1 _ g-2|a(t,nTo)P 

+ \l ^ |-)l/-(t)>, (26) 
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where 



where 



l±) 



V2 



(IT>±U», 



\f±{t)) 



^2(l±e-2|5(t,nro)P) 

•(|-(-l)"a(t,7iTo))±|(-l)"5(i,nro))). 



Vidal and Werner's analytic formula for negativity [17 
is used to find the negativity of such system, which is 



Vl - e-*|Q(t,nro)P 



(27) 



The negativity l|?7|) is plotted as shown in Fig. 
When the system just starts to evolve for some small 
t, there are some overlapping between the two coher- 
ent states associated with | t) and | i). As time 
goes on, the coherent states gets further apart and 
(— (— nro)|(— nro)]) « 0, and thus system 
become more entangled. The negativity approaches its 
saturated (maximum) value of i and the participation 
ratio 2. 

If the initial state is a thermal ground state at temper- 
ature T given in terms of density matrix as 



Pi{0) = 



IT) + U) (Tl + ai 

r=0 

T) + U) (Tl + (ir 



V2 V2 
(fa ■ Po{ct, ct*)\a){a\, 



(28) 



where PQ{a,a*) is the P-representation of the thermal 
state and is equal to 



1 _ Mf 

Po(a, a*) = ^^e 



(29) 



where the average boson distribution for the oscillator 
given by ([TS]). When T — 0, the state becomes the pure 
ground state |0)(0| which can be again manipulated as 
shown previously for a = 0. As in appendix \X[ by (jASp . 
the state can be calculated analytically at t = utq as 



p(nTo) 



T)(T 
T)(i 
i)(T 



(f'a ■ Pf|(a, a*, nTf))\a){a\ 
(fa ■ P-\i{a, a* , nTo)\a){a\ 
(fa ■ Pj^i (a, a* , nTo)\a) {a\ 
<fa ■ Pii{a,a*,nTa)\a){a\{ZQ) 



Pll{a,a*,nTo) 

Pl'^{a,a*,nTo) 
Pll{(^,a*,nTo) 



1 |Q + (-l)"2r.co|^ 

-e 



(31) 



2Trnr 

2'Knr 

e2(-i)"(irr+2)("-"*)"ao^ (32) 
2Trnr 

,g-2(-l)"(J;:+2)(a-Q-)nao^ (■33) 

(34) 



1 _ iQ_^^^iriii^ioi; 
-e 



2imr 



The participation ratios at i = titq with respect to the 
two-level system and the oscillator are then given respec- 
tively by 



Ka{nTo) 
Kr{nTo) 



1 + cxp [-16n2ao^(l + 2n^)] '' 
2(1 + 2nr) 



1 -I- exp 



(35) 
(36) 



As time goes on, K„ 2 and 2(1 + 2^^). The 

discrepancy comes from the fact that the state is a mixed 
state. 

The density matrix of the state can be computed at 
any time by solving the master equation (|13|) . This can 
be done numerically, but P|f(a,Q;*,i) and P||(a,a*,t) 
can also be computed analytically and so is K^, given by 



1 



:T(-l)("+^)5(t,nro)| 



21 



(37) 

where (5(t, nro) is defined in (|24p . The reduced density 
matrix of the quantum oscillator side is 

Pr{t) - I (fa-[P^^{a,a\t) + Pii{a,a\t)]\a){a\ 



(fa-P{a,a*,t)\a) {a\. 



(38) 



The participation ratio of the oscillator side can be eval- 
uated directly by 



1 



(fa / (fa'P{a,a*,t)P{a',a'*,t) 



— \a\'^ — \a'\^-\-a^a'+aa'* 



and is found to be 



1 -I- exp 



4|a(t,nTn)p 
2ftr + l 



(39) 



(40) 



It is plotted as shown in figure [TJ Putting t = utq in ((40|) 
restores (p6| . and putting T = (ft^ = 0) restores the 
zero temperature case in (j25p . With the numerical values 
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FIG. 1: Participation ratio of the oscillator Kr{t) of the sys- 
tem of the two-level system and the oscillator with coupling 
strength Ao = 0.2huio, with the ground state and the thermal 



state {j 



0.74239) as the initial states. 




10 to 



FIG. 2: Negativities J^{t) of the system of the two-level sys- 
tem and the oscillator with coupling strength Ao = 0.2fttJo, 
with the ground state and the thermal state (-^^ = 0.74239) 
as the initial states. 



of ps{t) computed, the negativity J^{t), a better measure 
of entanglement, is calculated using ([^0]) and plotted as 
shown in Fig. [2l It approaches the saturated value i as 
time goes on, but at a slower rate than the case for the 
coherent state as the initial state. 



V. ENTANGLEMENT WITH DECOHERENCE 

Quantum systems are vulnerable to decoherence due to 
its coupling with the thermal bath stated in ([4]), which 



is realistic since any quantum system cannot exist alone. 
To consider the the effect of decoherence, the state can 
be computed numerically by solving the master equations 
([T^ with all the terms. Its negativity can be computed 
numerically for every state by its definition ([20]). After 
the interaction with the bath, the entanglement of the 
system increases until some point and it will eventually 
decrease, as shown in Fig. |31 Although the negativ- 




FIG. 3: Negativities J^{t) of the system of the two-level sys- 
tem and the oscillator with coupling strength Ao = 0.2hijJo 



with the ground state and the thermal state (jp^r = 0.74239) 



as the initial states with decoherence parameters F 

0.01 

T"0 



c 



ity goes to zero at the end, it is still an open question 
whether the state is really non-entangled. Without de- 
coherence, the state approaches the maximum negativ- 
ity of ^ after a long time. However, with even a small 
amount of coupling, the state reaches a maximum nega- 
tivity which is less than i at some certain time and then 
damps eventually. 



The state affected by decoherence can be studied by 
finding the Glauber-Sudershan representation of the re- 
duced density matrix of the oscillator side analytically 
Its reduced density matrix can be found by trac- 
ing the spin components and be represented in Glauber- 
Sudershan P-representation as in ([55]) . Putting the ex- 
pression for the representation to the master equation 
the Fokker-Planck equations for the oscillator state 
are derived (in ignorant of the last term of (fT3|) ) 



5Pi 



TT/ii 



dt 



^—iujQt 



9Pi 



TT/ii 



da* 



da 



d 



d 



da da* 



d^P: 



TT/ii 



TT/ii 



dada* 



(41) 
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where the upper sign denote the term for and the 
lower sign for P||. The last term of (fT3)) just switches 
the spin at time t = nro where To = 5- Between these 
switches, the state can be fully described by the Fokker- 
Planck equations (|^T|) . The solutions to these two equa- 
tions can be found using the Green's function of the 
Fokker-Planck equation, as listed in the appendix [Bl If 
the initial state is a ground state, then the Glauber- 
Sudarshan representation of the oscillator at any time 
utq < t < (n + 1)tq is given by 



-PTT/ii("'"*'0 
1 



(42) 



27mj.(l 



exp 



.T(-l)"+'5(i,nTo)|2 



nr(l 



and if the initial state is a thermal state, then it is 



|aT(-l)"+^a(^,nTo)| 

fir 



(43) 

Here a{t,nTo) denotes value of the modal coherent state 
in any time t between htq and {n + 1)to, expressed as 



a{t, uTo) 



aoujQe 2 



(ri-l)To 



UJQ 



iC 
2 



{n(l 



_ giwo ( t - "To ) g - |- ro 



-(Mo + §)(t-nTo) 



(44) 



1 



Note that this a(i, nro) is different from another one 
for non-decoherence case in ([M)) . When C 0, 
a(nTo,7iTo) — » 2nao, which is the result found for the 
case without decoherence. However, for a non-zero C, 
a{t, tltq) for t oo. 

The representations in (|42p and (HS)) are important 
when studying the entanglement. By the partic- 

ipation ratio for the ground state is. 



Kr{t) = 



2 [2fv(l 



-ct\ 



1] 



1 + exp 

and that for the thermal state is 

Krit) 



4:\a(t.nTo) 1 2 
2n,.(l-e-c*)-|-l 



2 [277,^ -1- 1] 



1 -I- exp 



4|a(i,nTo) 



(45) 



(46) 



They are plotted as shown in Fig. ID At the beginning, 
the values of K are different since the ground state is 
pure at the beginning but the thermal state is already 
mixed. However, as time goes on, the two participation 
ratios coincides with each other. The reason is that the 
ground state is approaching the thermal state in order 
to achieve thermal equilibrium, as it can be seen in (I42|) 
that as < — *■ oo, it is asymptotically approaching (|43p . 

With ground state as the initial state, the maximum 
negativity as a function of the decoherence parameters 
(with F = C) is plotted as shown in figure [5l which 
shows that the entanglement decreases as the decoher- 
ence strength increases. 




FIG. 4: Participation ratio of the oscillator Kr{t) of the sys- 
tem of the two-level system and the oscillator with coupling 
strength Ao = 0.2fiti;o, with the ground state and the thermal 
state = 0.74239) as the initial states. The decoherence 

parameter is C — 0.1. 



A', 



-ToT 



FIG. 5: Maximum negativity Afmax of the system of the 
two-level system and the oscillator with coupling strength 
Ao = 0.2hujo with the coherent state as the initial states as a 
function of the decoherence parameters ( with F = C) . 



ToF 



FIG. 6: Time tmax (in terms of the unit of tq) to achieve max- 
imum negativity Mmax of the system of the two-level system 
and the oscillator with coupling strength Ao = Q.2huJo with 
the coherent state as the initial states as a function of the 
decoherence parameters (with T = C) . 
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I.Q3Q 

:,r=a.Q3i 
:,r=a.a3s 

■..T=a. 033 



t 



FIG. 7: Negativities J^{t) of the system of the two-level sys- 
tem and the oscillator with coupling strength Ao ~ 0.2fiti;o 
with the coherent state as the initial states with decoherence 
parameters T = C between Si^M and 




FIG. 8: A magnified graph of the negativities Af{t) of the 
system of the two-level system and the oscillator with coupling 
strength Ao = 0.2hu!o with the coherent state as the initial 
states with decoherence parameters F = C being equal to 
M31 and 



It is worth to note that the time i„ 



to achieve 



this maximum negativity decreases as the decoherence 
strength increases, with some jumps at a number of 
points as shown in figure [SI The jump can be explained 
by observing the time evolution of the entanglement with 
different values of r(= C). As F increases, the graph 
of the time evolution of the entanglement shifted down- 
wards continuously as shown in figure [T] Unlike the time 
evolution of the entanglement while there is no decoher- 
ence which increases monotonically (as in figure [2]), there 
are some ripples in the variation as time goes on when 
decoherence due to a thermal bath is present. The max- 
imum negativity Nmax still decreases continuously as F 
increases, but the time tmax to achieve this maximum get 
a jump because of these ripples, as shown in the magni- 
fied graph of the time variation of negativity in figure [S) 



VI. CONCLUSION 

It is shown that the Jaynes-Cummings interaction be- 
tween a quantum oscillator and a Pauli-like two-level sys- 
tem, with the help of a delta jump in the two-level sys- 
tem, enhances the entanglement between the two systems 
to its maximum possible values, no matter the system is 
initially a Pock state |0) or a thermal state. The coupling 
of each subsystem to a Markovian bath causes dissipation 
of the entanglement. The system will be in thermal equi- 
librium with the Markovian bath after a very long time. 
The maximum time for achieving the maximum entan- 
glement shows discontinuous jumps over the parameters 
of decoherence. 
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APPENDIX A: ANALYTIC EXPRESSION OF 
THE STATE WITHOUT DECOHERENCE AT 

t = nro 

Suppose the Dirac delta pulse ([3]) is imposed as the 
tunneling to the system and there is no decoherence. The 
time evolution between the pulses is given by 

C/i = D{a^a,)e-'^-''^D\aocr,), (Al) 

where D is the displacement operator given by 

L»(a) ^e""*-"*". (A2) 



The effect of the Dirac delta pulse on the system is merely 
{—iux), and thus the time evolution at time t = nro is 
given by 



U{nTo) = i-ia^UiT. 

By cr^D{aQ(7;,)a^ = D^aoa^) 

gi7rata^(Q,^^^-)g-™ata ^ D^aQa^), wc havc 



, , J D'^{2naoaz) lor even ti 

(nro) = I icr^e-'''''''''D^2naoa,) for odd n 



for even 



(A3) 
and 

(A4) 



As a result, if the system is a pure state desccribed by 
kets, the initial state \ipo) will evolve to 

U{nToM>). 

And if it is a mixed state, the state is described by density 
matrices and at t — utq the matrix is given by 



p(nro) 



D'' {2naoaz)poD{2naQaz), for even 



for odd n 



^'"Dt(2naotT.)po^(2nao<T.)e^™'v„ 



(A5) 
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APPENDIX B: GREEN'S FUNCTION OF THE 
FOKKER-PLANCK EQUATION 



P|l/^(a,a*,7iTo < t < (n+ 1)to) (B1) 
= J d^a ■ P{a,a* ,t\a' ,a'* ,nTo)P-\i{a' ,a'* ,nTo), 

where the Green's function is given by 

1 



P{a, a* , t\a , a *, htq) 



■ exp 



(B2) 



The function w{t2,ti) is actually given by the driving 
force 



dt' {^aoLL!Q)e 



iuJo{t2-ti-t') 



e 2 



iuJo{t2-ti)_ 



- {iujQ + §) 



for the upper sign denotes the spin-up and the lower sign 
the spin-down. Then it can be proved that 



The Fokker-Planck equations ((4T|) can be solved with 
the use of Green's function [19] . Suppose at time t = utq, 
the state is given by P|x (a, a* , nro), the state at any time where 
riTo <t< (n + 1)to is then given by 



w{t,nTo) = (-l)"e-^«^Oy,(i _ nTo,0), (B3) 



w{t — utq, 0) = ±(— l)"Q!o(j-'oe 



iujQt 



LOo - «§ 



(B4) 
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